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CN ■ Abstract 



As strong evidence for inflation, the relic gravitational waves (RGW) have been extensively studied. 
Although, it has not been detected, yet some constraints have been achieved by the observations. 
0> ' Future experiments for the RGW detection are mainly two kinds: the CMB experiments and the laser 

tN , interferometers. In this paper, we study these current constraints and the detective abilities of future 

experiments. We calculate the strength of RGW flg{k) in two methods: the analytic method and the 
numerical method by solving the inflationary flow equations. By the first method we obtain a bound 
QQ ' fig < 3.89 X 10~^^ at ^ = O.lHz, where we have used the current constraints on the scalar spectral 

lO ■ index, the tensor-scalar ratio, furthermore, we have taken into account of the redshift-suppression effect, 

^T ' the accelerating expansion effect, the neutrino damping effect on the RGW. But the analytic expression 

of fig(fc) depends on the specific inflationary models and only applies for the waves with very low 
frequencies. The numerical method is more precise for the high frequency waves and applies to any 



f^ . single-field inflationary model. It gives a bound fig < 8.62 x 10 ^^, which is independent of the 



inflationary parameters, and applies to any single-field slow-roll inflationary model. After considering 
f->i , the current constraints on the inflationary parameters, this bound reduces down to fig < 2 x 10^^^. 

These two methods give the consistent conclusions: The current constraints on the RGW from LIGO, 
big bang nucleosynthesis, and pulsar timing are too loose to give any stringent constraint for the 
single-field inflationary models, and the constraint from WMAP are relatively tighter. The future laser 



O 

j^ , interferometers are more effective for detecting the RGW with the smaller tensor-scalar ratio, but the 

CMB experiments are more effective for detecting the waves with the larger ratio. These detection 

methods are complementary to each other for the detections of RGW. 
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1. Introduction 

In the past, a number of observations on the CMB power spectra [1, 2, 3] and on the large scale structure (LSS) 
[4] have supported inflation as the good phenomenological model in the sense that it naturally gives rise to the 
origin of the primordial fluctuations with a nearly scale-invariant and gaussian spectrum. The overall expansion 
of the Universe at very early stage, as well as the evolution of fluctuations of the perturbed spacetime metric, can 
be accounted for in the framework of inflationary models. In addition to the primordial density perturbations, 
inflationary models also predict a stochastic background of RGW, which is the tensorial perturbations of spacetime 
metric. The detection of such a background would provide incontrovertible evidence that inflation actually occurred 
and would also set strong constraints on the dynamic of inflation [5] . 

There are mainly two kinds of experiments to detect RGW at correspondingly different frequencies. For RGW 
of very low frequencies v ^ 10^^^ ^ lO^^^Hz, one can observe them by detecting the power spectrum of CMB 
B-polarizations [6] . Now, the first-three- year results of WMAP [2] have not yet found the evidence of RGW. The 
experiment of the Planck satellite [7] with higher sensitivity to polarizations is scheduled for launch in 2007, and 
the Clover (Cl-Observer) [8] and CMBPol [9] projects with much higher sensitivities than Planck arc also under 
development. For RGW of high frequencies v ~ 10^'' ^ lO^Hz, another kind of experiments apply, i.e. the laser 
interferometers detectors, including the current TAMA [10], VIRGO [11], LIGO [12, 13], and the future LISA [14], 
ASTROD [15], BBO [16], and DECIGO [17]. Besides these two kinds, other methods have also been used to constrain 
the strength of RGW. For example the timing studies on the millisecond pulsars, which can constrain the amplitude 
of the gravitational waves by studying the signal residuals of the millisecond pulses [18]. This method is sensitive 
to the waves with frequencies at 10~^ ~ 10~^Hz. The electromagnetic detectors are based on the processes of 
resonant responses of electro-magnetic field to incident gravitational waves [19, 20, 21]. This kind of detectors are 
aiming at detecting RGW at very high frequencies v ~ 10^ — 10^°Hz, which can only be produced by inflation [22]. 
The observational results of BBN also can constrain the strength of RGW [23, 24, 25] at all frequencies. Although 
the RGW have not been found yet, some constraints on it have already been obtained by these experiments or 
observations. 

This paper is to study the various modifications on the power spectrum of RGW and examine the constraints on 
it from the experiments. On the spectrum, we will consider the modifications due to such important effects as the 
redshift-suppression by various periods of cosmic expansions, the current accelerating expansion, and the damping 
by free-streaming of neutrinos. Two methods will be applied in this study: the analytic method and the numerical 
method by solving the inflationary flow equations. After considering all these damping effects, we will arrive at an 
analytic formula of the strength of RGW, which, as a function of the wavenumber fc, depends on the scalar spectral 
index n^ and the tensor-scalar ratio r. By taking into account of the current observational constraints on Ug and r, 
we will obtain an upper limit of the strength of RGW fJg < 3.89 x 10^^^ at z/ = O.lHz. From the plots in the r — fJg 
plane, we find that the BBO experiments can detect the RGW if r > 8.3 x 10"'^ is satisfied. This means that BBO 
is more sensitive than Planck satellite, but less than Clover and CMBPol. But the would-be ultimate DECIGO 
can even detect for r > 6.8 x 10^®, much more sensitive than the CMB experiments. In this analytic method, fig 
depends explicitly on the ratio r, an undetermined parameter, whose value varies for the various specific inflationary 
models. Furthermore, the approximate power law of primordial power spectrum also may yield fairly large error if 
it is extended over a larger range of frequencies. 

To overcome these shortcomings, we move on to the numerical method. Because the RGW depends more 
sensitively on inflationary stage, during which it is generated, and there are a number of inflationary models, to take 
care of the predictions from these models, the numerical method is used, by which a great many of realizations are 
produced, representing the respective inflationary models. The inflationary flow equations are applied to numerically 
calculate the RGW, whereby an upper limit ilg < 8.62 x 10^^'* is obtained for any slow-roll single-scalar-field 
inflationary model, independent of any inflationary parameters. By taking into account of the current observed 
constraints on n^, a, and r, a much tighter limit Og < 2 x 10~^^ is arrived at r ~ 0.03, which is beyond the sensitive 
range of BBO. In numerically generating 10^ realizations, we flnd all of them satisfy the current constraints on ilg 
from LIGO, from pulsar timing, and from BBN, but only nearly 0.05% of them satisfy the current constraints on n,,, 
a and r. From the resulting r — f2g plane, one flnds the DECIGO, if put into running, will be effective for detecting 
the RGW with smaller r, but the CMB experiments, such as Planck, Clover and CMBPol, are more effective for 
detecting the RGW with larger r. They are complementary to each other for RGW detections. Our result from 
this numerical investigation applies only to the single-field inflationary models with the chosen initial conditions of 
Hubble slow-roll parameters. 

The organization of this paper is as follows. Section 2 gives a simply review on the RGW and its evolution 
equation. In section 3, an analytic expression of the strength of RGW will be obtained with the three damping 
factors being included, presenting the modiflcations due to the mentioned effects. Section 4 is devoted to the 
numerical computations, where the strength of RGW is numerically calculated by solving the inflationary flow 



equations. Finally section 5 is the conclusion. 

2. The relic gravitational waves and their evolutive equation 
Incorporating the perturbations to the spatially flat Robertson- Walker (FRW) spacetime, the metric is 

ds^ = a{Tf[dT^ - ((5y + h,j)dx'dx3] , (1) 

where a is the scale factor of the universe, r is the conformal time, which is related to the cosmic time by adr = dt. 
The perturbation of spacetime hij is a 3 x 3 symmetric matrix. The gravitational wave field is the tensorial portion 
of hij, which is transverse-traceless dih^^ — 0, S'^^hij = 0. Since RGW is very weak, \hij\ ^ 1, one needs just study 
the linearized evolutive equation: 

d^iV^d^h,,) = 16^Ga2(T)n,, , (2) 

where Hij is the tensor part of the anisotropy stress, satisfying 11^^ = 0, and 9^11^ = 0. It couples to hij as an 
external source. In the cosmic setting, 11^^ can be generated by the free-streaming relativistic particles [26, 27], the 
cosmic magnetic [28], etc. It is convenient to Fourier transform these quantities as follows 

h,,ir,^)=y^Vm^ f -^e['\k)h^{r)e^- , (3) 



(27r)3/2 '^ 



n,,(r,x) ^ V^/l6^ / -A^e\'\k)Ilt{T)e^^- , (4) 



where the index A = "+" or " x " labels the two polarization states of the gravitational waves. The polarization tensors 
e], are symmetry, transverse-traceless k^e\, (k) = 0, (5*^e^- (k) = 0, and satisfy the conditions e^^'"^^ {'k)e-- (k) = 
2S\\' and e^- (— k) = e^-- (k). Since the RGW is assumed to be isotropy and each polarization state has the same 

evolution and gives the same contributions, h\^ (r) is denoted simply by hkir), and U^ (r) by nfe(r), where fc = |k| 
is the wavenumber related to the frequency by i^ = k/27T (the present scale factor is set ao = 1). Then Eq.(2) can 
be rewritten as 

hk + 2-hk + k^hk = 16TrGa^{T)Ilk (r) , (5) 

a 

where the overdot denotes a conformal time derivative d/dr. Usually the interactions between gravitational waves 
and other matters are very weak, in many cases, the source Hj, in Eq.(5) is negligible, and the evolution of RGW 
only depends on the scale factor and its time derivative. But in this paper we include this source term, so that the 
damping from neutrino free-streaming is properly taken care of. 



3. The analytic power spectrum of RGW 

The primordial power spectrum of RGW 

As mentioned in the introduction, inflationary expansion, as an attractive idea to describe the very early universe, has 
received strong support from the observations of CMB anisotropies and from studies of the large-scale distribution 
of galaxy. There have been a number of models proposed. Here we will consider only the single field models. In the 
context of slow-roll inflationary models, the observables depend on three slow-roll parameters [29] 

.„.M(;)\ ,.,.„,(-), ,..M,(i^), ,a, 

where Mpi = (SttG)^^" — inp\/y8TT is the reduced Planck energy, V{(p) is the inflationary potential, and the prime 
denotes derivatives with respect to the field cj). Here, ey quantifies "steepness" of the slope of the potential, rjv 
measures "curvature" of the potential, and ^y quantifies the "jerk" . All these three parameters must be smaller 
than 1 for inflation to occur. One of the important predictions of inflationary models is the primordial scalar 
perturbation power spectrum, which is nearly gaussian and nearly scale-invariant. This spectrum can be written in 
the form 



Psik) = Psiko) { ^J , (7) 



where Hs is the scalar spectral index, a = dus/dlnk is its running, and fco is some pivot wavenumber. In this 
paper, ka = O.OSMpc"^^ is taken, which corresponds to a frequency ~ lO^^^Hz. The observations of WMAP gives 
Ps{ko) ~ 2.95 X 10^^A{ko) and A{ko) = 0.9±0.1 [1]. Another major prediction of inflationary models is the existence 
of RGW. The primordial power spectrum of RGW is defined by 

PT{k) ^ ^^^^h+hk , (8) 

TT 

where hk is the solution of the equation in (5). This spectrum can also be put in a simple form 

J, \ nt(feo) + ^at ln(fc/fco) 



PTik)^PT{ko)l^-j , (9) 

where nt{k) is the tensor spectral index, and at = dnt/dlnk is its running. In the single-field inflationary models, 
a standard slow-roll analysis gives the following relations 



nt 



n 8 16 

K-l) + gJ , r ^ -{I - n,) + —7JV , (10) 



where r{k) = PT{k) / Ps{k) is the so-called tensor-scalar ratio. These formulae relate the tensorial parameters rit and 
at to the scalar parameters n^ and r; the latter are accessible to the observations of CMB and of LSS. As shown in 
Eq.(lO), the relation between r and Ug involves the slow-roll parameter rjv, depending on the specific inflationary 
potential. Inserting these into Eq.(9), one has 

PT{k) ^ Ps{ko) y^ r X f — j . (11) 

In general, the tensor-scalar ratio r may vary with the wavenumber k. Here and in the following sections we will 
take the value of r at fc = fcoj i-e. r = r(fco). Now the primordial spectrum of RGW only depends on the parameters 
Us and r. The recent constraints by the observations of three-year WMAP, SDSS, SNIa and galaxy clustering [30] 
are 

Us = 0.965 ± 0.012 , (68% C.L.) (12) 

r < 0.22 , (95% C.L.) . (13) 

From the derivation it is obvious that the formula (11) applies properly for in a frequency range around ^ lO^^^Hz. 
The strength of the gravitational waves can be also characterized by the gravitational waves energy spectrum 

where pc = 3Hq/8ttG is the critical density and Hq = lOO/i km s~^Mpc~^ is the present Hubble constant (the value 
h = 0.72 is taken through this paper), fig can be related to the primordial power spectrum by the formula [27, 31] 



1 

^0 



^aik) = j^j^k^PT{k)T^{k) , (15) 



where the transfer function T(k) will take into account of the various damping effects mentioned early, and will be 
discussed below. 

The damping effects 

Here three kinds of damping effects will be addressed: First we only consider the redshift-suppression effect caused 
by the overall expansions of the universe. So temporarily we drop the the anisotropy stress term Hfe(T) in Eq.(5) 
due to neutrino free-streaming, 

hk + 2-hk + k^hk = . (16) 

a 

This equation of RGW only depends on the behavior of the scale factor air). It has been known that, during the 
expansion of the universe, the mode function /ife(r) of the gravitational waves behaves differently in two regimes [22]: 
far outside the horizon {k <C aH), and deep inside the horizon (fc 3> aH). When waves are far outside the horizon, 
the amplitude of hk keeps constant, and when inside the horizon, the amplitude is damping with the expansion of 
the universe 

hk « —TT ■ (17) 



By numerically integrating the Eq.(16), this effect can be approximately described by a transfer function [32] 



where keq — O.OTSOm/i^Mpc^^ is the wavenumber corresponding to the Hubble radius at the time that matter and 
radiation have equal energy densities. And ro = 1.41 x 10''Mpc is the present conformal time. It is obvious that, 
this factor ti{k) is oscillating with wavenumber k due to the Besscl function ji(fcTo). In practice, one is usually 
interested in the overall outline of the amplitude of RGW, as the quick oscillations with k are of no importance. For 
the waves with fcro 3> 1, this factor can be written as 



The above transfer function (19) does not include the effect of accelerating expansion of the present universe, 
which has been indicated by the observations on SN la. The spectrum of RGW has been studied in specific 
models for dark energy [33], such as the Chaplyngin gas models and the X- fluid model. In the Ref. [22], we have 
presented an analytic solution of RGW in ACDM universe, and found that the amplitude of the gravitational waves 
has been modified by the presence of the dark energy during the current expansion. In the higher frequency range 
(i' ^ 3 X 10~^^IIz) that we are interested in this paper, the amplitude acquires an overall factor Q.m/^h as compared 
with the decelerating model, where f2m and IIa are the present energy densities of matter and vacuum, respectively. 
So this effect can be simply described by a damping factor, 

'- §7 ■ (20) 

In the standard ACDM model with Vim = 0.27 and 51a = 0.73, this effect contributes a damping factor of t^ ^ 0.137 
for the strength of RGW in Eq.(15). 

The third to be considered is the damping effect of the free-streaming neutrinos [26], i.e. the anisotropic stress 11^ 
on the right-hand of the Eq.(5). This effect has been considered by Weinberg [26]. This effect is primarily produced 
by neutrinos during the radiation dominated era when they are decoupled and are free streaming in the universe, 
especially right after the gravitational waves enter the the horizon. The overall amplitude of RGW will reduce roughly 
by an amount of 20%. It has been shown that anisotropy stress can reduce the amplitude for the wavelengths that 
re-enter the horizon during the radiation-dominated stage, and the damping factor is only dependent on the fraction 
/ of the free-streaming relativistic particles over the background (critical) energy density in the universe. For the 
waves that enter the horizon at later times, the damping effect is less important. A number of works have done to 
discuss this effect, and Reference [27] has found that the effect can be approximately described by a transfer function 
ts for the waves with wavenumbers k > lO^^^Hz (which re-enter the horizon at the radiation-dominant stage), 

_ 15(14406/4 - 55770/3 ^ 3152975/2 _ 48118000/ + 324135000) 
*^ ~ 343(15 + 4/) (50 + 4/) (105 + 4/) (180 + 4/) ' ^ ' 

When the wave modes (10~^^Hz< k < 10~^°Hz) re-enter the horizon, the temperature in the universe is relatively 

low (< IMeV), the neutrino is the only free-streaming relativistic particle. Thus we choose the fraction / = 0.4052, 

corresponding to 3 species of neutrinos in the standard model, the damping factor is then tg, — 0.80313. But for 

the waves with very high frequency (fc > 10^ ^"Hz), the temperature of the universe is still very high when they 

re-enter the horizon, the neutrinos were still in interaction, and the value of / is quite uncertain. This is because 

the detail of how many species of particles are free is not accurately known and depends on the cosmic environment 

and the particle-interaction models. Thus, the detection of RGW at this range of frequencies offers the possibility 

of learning more about the free-streaming fraction / in the very early universe. So we choose / = 0, i.e. there are 

no free-streaming relativistic particle at that time. The corresponding factor is then ^3 =: 1. But for the waves with 

k < 10~^^Hz, which re-enter the horizon during the matter-dominated stage, the neutrino density is so small that 

its damping impact can be neglected, so we choose t^ = 1. Overall, the neutrino damping effect can approximately 

be summarized as 

r 1, fc < IQ-^'^Hz 

<3(fc) ~ <^ 0.80313, 10-^^Hz <k< 10-^°Hz (22) 

[1, fc > 10-^°Hz. 

Putting these three effects together, the total transfer function is the combination of three factors 

T{k) = ti X t2 X ta , (23) 



among these ti is dominantly important, which approximately shows the evolution of RGW in the expanding universe. 
The function ^2 has the relatively smaller damping on RGW, which accounts for the accelerating expansion of the 
universe quite recently {z ~ 0.3). The value of energy spectrum f2g is reduced by nearly an order by this effect. 
The function of is has the most uncertainty, as has been discussed above. In the extreme case with / = 0, one has 
ts — 1, i.e. no damping; and in another extreme case with / = 1, one has ts = 0.35, the smallest value of ^3. In the 
case oi f — 0.4052, ^3 = 0.80313, only contributing a damping factor t| = 0.645 for the strength of the RGW. 

There are some other possible mechanisms, which might affect the amplitude of the gravitational waves. For 
example, the QCD transition [34, 35], e+e~ annihilation [34, 36, 35], the cosmic reheating [37, 22], and so on [27]. 
These could influence the value of expansion rate a/a, and therefore affect the strength of RGW. However, these 
effects are either small, as shown in literature, in comparison with the effects we have discussed, or there are some 
uncertainties in their analysis, so these are not considered in here. 

The upper limit of fig and the sensitivities of future experiments 

The future detectors of RGW arc mainly classified into two kinds: one kind is through CMB for very low 
frequencies, and another is based on laser interferometers for relatively high frequencies. 

For the waves of very low frequencies i^ < 10^^'"' Hz, the CMB experiments are sensitive. For instance, the Planck 
satellite can detect the RGW if r > 0.1 [7], the ground-based experiment Clover can detect the signal if r > 0.005 
[8], and CMBPol can detect if r > 10~^ is satisfied [9]. It should notice that if r < 1 x 10"**, the RGW may not be 
detected by the CMB experiments. This is because the CMB B-polarizations generated by the cosmic lensing are 
also very large, and the signals from the RGW may be subdominant to the lensing effects. [38]. 

The direct detections the RGW by laser interferometers are sensitive to the waves with high frequencies. For the 
waves with k > 10^*"Hz, inserting the formulas (19)-(23) with ^3 = 1 in Eq.(15), the strength of the gravitational 
waves becomes 

22.5 Prik) fn. 

Using the expression of Prik) in Eq.(ll), one gets 



"^^'^^I^#i;lTf) - 1-08 X 10-P.W , (24) 



-i+f6[("=-i)+i]in(fe/fco) 
ng{k) c^ 2.87 X 10-^^ r { y ) , (25) 



k 



where ^(fco) — 0.9 has been taken. This function depends on the wavcnumbcr fc, the tensor-scalar ratio r and the 
scalar spectral index n^. It should be pointed out that, just as Eq.(ll), this formula (25) also applies properly in a 
frequency range around ~ lO^^^Hz. 

The advanced LIGO can detect the waves with flgh^ > 10^^ at :^ ~ lOOHz [13]. The LISA project is expected 
to detect waves with ilgh'^ > 10^*^ at i^ ~ 0.005Hz [14]. The ASTROD, a space project sensitive to the waves with 
frequency at i^ <E (10^^, 10^^)Hz [15], is expected to detect the waves with flgh"^ > 10^*^ at t^ ~ 5 x lO^^Hz. The 
BBO, another important project, can detect a background RGW with fig > 2.2 x 10^^^ at i^ ~ (O.l-l)Hz [16]. The 
DECIGO project, having a much higher sensitivity by design, is expected to detect the RGW with ftgh^ > 10^^° 
at i^-O.lHz [17]. 

First, we will estimate the upper limit on the strength of RGW in Eq.(25). Here we assume Ug < 1 and r < 0.22, 
which are consistent with the current observations [30]. The formula (25) gives an upper limit of fig at j/ = O.lHz: 

Qg < 3.89 X 10-1*^ . (26) 

And this limit is arrived at n^ = 1 and r = 0.22. This limit is nearly an order smaller than the result in Ref.[31]. This 
is because our analysis has taken into account of the damping effect of the accelerating expansion of the universe 
and the running of nt in the primordial spectrum. This limit is in the sensitive ranges of BBO and DECIGO, but 
beyond those of LIGO, LISA and ASTROD. In Fig.[l], we plot the strength of RGW at 1/ == O.lHz, as the function of 
r, where several models with different n,, are demonstrated. One sees that when r < 0.01, the curves of the function 
fig are almost overlapped for the models of different rig, and only depend on the variable r. But when r > 0.01, the 
models of different n^ can be distinguished. For a fixed r, a larger rig yields a larger Qg. This figure also tells that 
BBO can detect the RGW if r > 8.3 x 10^^, so it is more sensitive than the Planck satellite, but less than Clover 
and CMBPol. It is interesting to notice that DECIGO can detect the RGW if r > 6.8 x 10~^, which is much more 
sensitive than all the CMB experiments (r > 10"^). 

The predictions of inflationary models 

The strength of RGW in Eq.(25) depends on the values of n^ and r. Observations have yielded quite solid 
constraints on Ug, but the value of r is still uncertain. The relation between Ug and r depends on the specific 



inflationary models, and diflFcrcnt models will predict very different r. In the following examinations will made for 
several inflationary models, which predict different values of r. One may categorize slow-roll models into several 
classes according to the parameter space spanned by n^, a and r [39]. Each class should correspond to specific 
physical models of inflation. Here we categorize the models according to the curvature of potential rjv in Eq.(6), as 
it is the only parameter that enters into the relation (10) between ng and r. The classes are defined in the following: 

Case A: negative curvature models rjy < 

The negative ryy models often arise from a potential of spontaneous symmetry breaking. One type of often- 
discussed potentials have the form of V = A'' [1 — (0//i)P], where p > 2. This kind of models predict the red tilt 
rig < 1, which is consistent with the observations of three-year WMAP. Also these models predict pretty small r. 
For the model with p = 2 in Ref . [39] , 

r~8(l-n,)e-^(i-"=) , (27) 

where N is the number of e-folds, taken be in the range N G [40, 70] to account for the current observations on 
CMB [40, 1, 2]. Here we choose the value N = 70. Using the constraint on rig in Eq.(12) yields the constraint 
r e [0.014, 0.037]. From Fig.[l], one finds this is beyond the sensitive range of the Planck satellite, but in the 
sensitive ranges of Clover, and CMBPol. And it is also in the sensitive range of BBO and DICIGO. In other models 
with p > 2, the predicted values of r arc much smaller than that of the model with p = 2. 

Case B: small positive curvature m,odels Q < rjv < 2ey 

These models contain as two examples the monomial potentials V — h'^{(j)/ ii)^ with p > 2 for < 77^ < 2ev and 
the exponential potential V = A"^ exp{(j)/n) for rjv = 2ey. In these models, to the first order in slow roll, the scalar 
index is always red rig < 1 and the following constraint on r is satisfied 

^(l-n,)<r<8(l-n,) . (28) 

Using the constraint on n,, in Eq.(12), one finds that r E [0.061, 0.376], which is in the sensitive ranges of Clover, 
CMBPol, BBO, and DECIGO. The sensitivity limit of Planck is just in this span, so it may be able to detect the 
model. 

Case C: intermediate positive curvature models 2ey < rjv < 3ey 

The supergravity-motivated hybrid models have a potential of the form V ~ A^ [l + a\\i{cf)/Q) + A(0//i)''] , up 
to one-loop correction during inflation. In this case, 

n, < 1 , r > 8(1 - n,) , (29) 

are satisfied. Using the constraint on n., in Eq.(12), one finds that r > 0.184, which is very close to the current 
upper limit r < 0.22. Fig.[l] shows that this model is in the sensitive range of Planck satellite. 

Case D: large positive curvature models rjv > 3ey 

This class of models have a typical monomial potential similar to the Case A, but with a plus sign for the 
term {(J>/ij.)P: V^ = A^ [1 + {(()/ ii)'p]. This will enable inflation to occur for a small value of < mpi. This model 
predicts a blue tilt of scalar index n^ > 1, which is contradict to the constraint in Eq.(12). But we should notice 
that the observations of three-year WMAP has not yet ruled out the blue spectrum. If the the running of Ug with 
the wavenumber k is allowed, then the best fit of WMAP data suggests that ns{k = 0.002Mpc^^) = 1.21^0 :[g, and 
a(fc — 0.002Mpc~^) = — 0.102^gQ|3[2]. This is a blue spectrum with a negative running. So the determination of 
the value of Ug depends on the more precise observations. 



4. The inflationary flow equations and the predictions for RGW 

In the discussions, the RGW given by the analytic expression (25) depends on the value of tensor-scalar ratio r, which 
has not yet been determined by the observations. Moreover, Eq.(25) is a good approximation only for the waves with 
wavenumber around fc ~ ka(^ 10~^^Hz). Therefore, for the waves of high frequencies, say with v ~ O.lHz, nearly 15 
orders larger than the value of fco, direction application of the formula (25) may lead to large errors. Consequently, 
it will be restricted in practice. To avoid these shortcomings, in this section, we will employ the technique of the 
inflationary flow equations to relate RGW in lower frequencies to that in higher frequencies. 

The inflationary flow equations 

The inflationary flow equations were first introduced by Hoffman and Turner [41] as a way of generating a large 
number of slow-roll inflationary models to be compared to the observational data. This method applies to any 
slow-roll single scalar field inflationary models, and relies on defining a set of Hubble slow-roll parameters, which are 



the derivatives of the Hubble parameter during inflation. The major advantage of this method is that it removes the 
field from the dynamics, and allows one to study the generic behavior of slow-roll inflation without making detailed 
assumptions about the underlying particle physics. In this section, we will also use this method to generate a large 
number of inflationary models, the observables of which are required to be consistent with the current observational 
constraints in low frequencies. Then we will numerically solve the strength of RGW in very high frequencies. In this 
method the Hubble slow-roll parameters are deflned by 
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where primes are derivatives with respect to the scalar field (p, and H{(j}) is the Hubble parameter as the function of 
(/), related to the potential V{4)) by the so-called Hamilton- Jacobi formula. 
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These Hubble slow-roll parameters satisfy an infinite set of hierarchical differential equations, called inflationary flow 
equations: 
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where TV is the number of e-folds of the inflation, and a = 2Ai - 
these flow equations: One is that e = 0, A; = for / > 2, and a~ 
only if cr > 0, this fixed point is stable, i.e. the attractor solution. 



4e. There are two families of fixed points of 
constant. In Ref.[42], the authors found that, 
The other family of fixed points arc given by: 



e=constant, a = — 2e, A2 = e , and A; = eA(_i for / > 3. Later we will show that the second family of fixed points 
is not stable. The slow-roll parameters tend to run to the attractor with the expansion of the universe, as long as 
the slow-roll condition e < 1 is satisfied. In order to actually solve this infinite series of equations numerically, it 
must be truncated at some I by setting a sufficiently high slow- roll parameter to zero, i.e. Am+i=b, with b being a 
constant, and Xm+2=0 for some suitably large m. In this section, we make the truncation of this series at m=10, 
and choose a set of acceptable initial conditions as in Refs.[42, 31]: 



XiU e 



eUe [0,0.8] , 

ct|, e [-0.5,0.5] , 
A2I, e [-0.05,0.05] , 
-0.025 X 5"'+^, 0.025 x 5"'+^] , (3 < I < 10) , 
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where the subscript \i denote the corresponding initial values. This set of eleven equations in Eqs.(32)-(34) is an 
autonomous system [42]. We choose the constant b ^0, and set the left hand of these equations to be zero. Then 
the only real solution for this eleven-equation set is given by 



&'/" , ^c = -2fei/" , A,, = 6'/", (2 < ? < 10) , 



(39) 



where the subscript c means the fixed point. This is just the second family of fixed points with Ec — b^^^^. As usual, 
in order to study the stability of this fixed points, let us consider the small perturbations around the fixed point, i.e. 

e = ec + Se, a = (7c + S<T, Xi= Xic + SXi, (2 < I < 10) . (40) 

Substituting these into Eqs.(32)-(34), one gets the first-order differential equations for the small perturbations 



(41) 





/ Se \ 




/ 5e \ 


d 

dN 


6a 




5a 




= M 






\ SXio J 




\ 5Xw J 



where the matrix M depends upon the values of ec, Uc and A/c, I = 2, ..., 10. If this fixed point is stable, at least, it is 
necessary that the real parts of the eigenvalues of the matrix M are negative [43] . We have performed calculations 
and found that no matter what value of 5 7^ is chosen, this condition can not be satisfied. Therefore, this fixed 
point is not stable. If we choose the value 6 = 0, this eleven-equation set in Eqs.(32)-(34) has fixed points, which 
belong to the first family mentioned early and are stable only if CTc > is satisfied. 

It is obvious that the evolutions of this eleven-equation set will be different for the conditions with different b. 
Although whether the fixed points are stable or not depends on the value of 6 = 0, still in the computation below, 
we will do computations for both kinds of initial conditions: one with 6 = 0, and the other 6 7^ with 



6e [-0.025 X 5"^ 0.025 X 5"*^] 



(42) 



It will turn out that, for these two cases, our calculational results of fJg are very similar. So in the following sections, 
we will only present the results of the 6 = case. 



The inflationary parameters and the strength of RGW 

Many observable parameters in inflationary models can be related to the . Here we are only interested in three 
such kind of observable parameters for the slow-roll inflationary models: the tensor-scalar ratio r, the scalar spectral 
index n^, and its running a. They are related to the Hubble slow- roll parameters as the following (up to the second 
order in the slow-roll) [44] 
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1 



(5- 



- 16e[l 
3c)e2 - 



1 dn. 



5c)(T€ 



> 



C)X2 



(44) 

where c = 4(ln2 + 7) — 5 ~ 0.0814514 (with 7 the Euler-Maschcroni constant) is a constant. Once the inflationary 
fiow equations in (32)-(34) are numerically solved, the values of these three observables are obtained. As we have 
pointed out early, we are interested for the gravitational waves in a very wide frequency range, i^ e [10~^^, 10^]Hz, 
and the primordial power spectrum in Eq.(9), as an analytic approximation, may not apply properly. Therefore, we 
need to adopt the following primordial power spectrum [45] 
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which is proper for the slow-roll inflationary models in general. Here H is the Hubble parameter of inflation when 
the waves exactly crossed the horizon with k — aH. If one ignored the small slow-roll parameter e in Eq.(46), then 

, a result for the exact de Sitter inflation, depending only on the Hubble 
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one would end up with Prik) = — 

parameter H. The formula (46) can be rewritten as 
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where e^ and Hi are the respective values of e and H when fco just crosses the horizon at a = k^/ Hi. By the way, 
the RGW power spectrum can be related to the scalar one by Prika) = Ps{kQ)r{kQ) as before. The value of H is 
also related to Hi through the parameter e by the following 



H{N) = H, exp 
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where Ni is the number of c-folds if iJ = Hi. Inserting the Eqs.(47) and (48) into Eq.(15), one obtains the energy 
spectrum of RGW 
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where T{k) is the damping factor, and Hq is the present Hubble constant. Making use of T{k) as given in Eqs.(19), 
(20), (21) and the ratio r in Eq.(43), one converts the above energy spectrum of RGW into the following form 
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which now depends only on the slow-roll parameters e and a. As an advantage to Eq.(25), this spectrum is good 
also for the waves with v ^ lO^^^Hz. Before numerically computing it through the inflationary flow equations, we 
first give an estimate of its upper limit. Since < e < 1 is satisfied during the inflation, Eq.(50) yields the upper 
limit when taking e = 

^Ak) < 7.34 X io-i3!l^£(^i!l±^ , (51) 

depending on the values of Ci and ai. When e^ = 1, the right-hand of this inequality has the maximum value, so one 
can give a loose upper limit of i^g{k): 

rigik) < 8.62 X lO"^"* , (52) 

where the approximation r ~ 16e has been used, and the second-order terms of the tensor-scalar ratio r have been 
omitted. This upper limit holds only if the slow- roll condition is satisfied. Compared with the limit in Eq.(26), this 
upper limit is arrived at without explicitly using the values of n,, and r. Besides, this limit applies for a wider range 
of frequencies i^ > 10~^"Hz. It is seen that this limit is much less stringent than that in Eq.(26). Still the limit (52) 
is much smaller than the sensitivity of, and therefore can not be directly detected by LIGO and LISA, however, it is 
now larger than the sensitivity of, and can be detected by ASTROD, BBO and DECIGO. In writing down the limit 
(52), ^3 = 1 has been used, which is valid for waves with i/ ^ 10~^''Hz. But for the waves with v S [10~^^, 10~^'']Hz, 
one should use t^ ~ 0.80313, and thus the limit becomes rig(fc) < 5.56 x 10~^^, which is a little tighter than the 
limit in Eq.(52). 

The current constraints on the cosmic parameters 

Through the discussion above, we know that the values of the inflationary parameters Ug, a, r, and fig are all 
directly related to the Hubble slow-roll parameters. Here, we give a review of the current constraints on them. The 
constraints on rig , a and r mainly come from the observations on large scales, including the observations of CMB, LSS 
and so on. Here we call them as "large-scale constraints" (LSC). When solving the inflationary flow equations, we will 
take the initial condition at the time of horizon-leaving, i.e. fco = aH, where the povit wavenumbcr feg = 0.05Mpc~"'^ 
as before. 

Now the WMAP CMB data (1st year) gives [1] n^ = 0.93 ± 0.07, a = -0.047 ± 0.040, and the best fit of 
WMAPext+2dFGRS galaxy survey gives n^ = 0.93 ±0.03, a = -0.03lto:oi7. A fit using WMAP CMB data and the 
SDSS galaxy survey gives [46] Ug = 0.98±0.02, a = -0.003±0.010. Combining the observations of three-year WMAP, 
SDSS, SNIa and galaxy clustering [30], one can give the constraints Ug ~ 0.965 ±0.012, a = — (2.0± 1.2) x 10^^ and 
r < 0.22. These bounds of Ug and a are all at the 68% confidence level, and that of r is at 95% confidence level. 
These various bounds are consistent with each other, taking into account of the corresponding confidence levels. 
Here in our calculation we will choose the most loose constraints 

Ug e [0.86, 1.00] , a e [-0.087, 0.007] , r < 0.22 . (53) 

which imply that the primordial scalar spectrum is "red" or scale-invariant, and the running of scalar index is very 
small, as required by the slow-roll infiationary models. 

The constraints on Qg mainly come from the observations at small scale. Here we call them as "small-scale 
constraints" (SSC), which include the tightest constraint from the observations of the pulsar timing [18] 

flgh^ < 2 X lO"'', ^ = 1.9 X IQ-^Hz ; (54) 

the constraint from the recent observations of LIGO [12], 

^g < 8.4 X 10"^ 69Hz <v < 156Hz; (55) 

and the constraint from the observations of BBN [23, 24], 

Wgh^ < 8.9 X 10"^ , (56) 

where Wg = J ng{h')dlni>. Comparing with the constraints (26) and (52), it is fair to say that the current SSC 
are too loose to give any constraint on the single-field inflationary models. This result will also be checked in the 
following numerical calculation. 

The distribution of the realizations 

In this part, we will present a numerical program to solve the inflationary flow equations (32)- (34) for a large 
number of models, where each initial condition randomly chosen within the constraints of (35)- (38) will represent 
an inflationary model. 



First, we want to study how tight the LSC of Eq.(53) and SSC of Eqs.(54)-(56) are, as constraints, on the 
inflationary models. We have produced 10^ reahzations of inflationary models. It turns out that all these realizations 
satisfy the SSC, which attests to the conclusion before: the current SSC are too loose to give any actual constraint on 
the single-field inflationary models. On the other hand, among these 10^ realizations, only 5523 of them (^ 0.05%) 
satisfy the LSC of Eq.(53). So this constraint is tighter for the inflationary models. In the following, we will mainly 
discuss the distribution of these 5523 realizations. 

During the numerical calculations, the inflation can end in one of the following two ways. One is that e < 1 is 
violated in the process of computing, then the inflation automatically stopped. A number of of inflationary models 
are of this class, such as the polynomial "large-field" models, the "small-field" polynomial potentials [39]. The other 
way is by an abrupt termination, perhaps from intervention of an auxiliary field as in hybrid inflation. The linear 
potentials and the exponential potentials also belong to this class [39]. Here we choose the abrupt stop to be at 
iV = 70 in computation. We have found that among these 5523 realizations, only 14 of them stop the inflation in 
the first way, all others do in the second way. This fact is consistent with the previous works [42, 31]. In following 
we will discuss these two kinds of realizations separately. 

First, we discuss the 5509 realizations that exit inflation by abrupt termination. They also satisfy both the large 
and small scale constraints in Eqs.(53)-(56). These models and can inflate at least 70 e-folds. In Fig. [2], we plot 
them in the r — fig plane, which shows the following characters: 

a. For a fixed r, the distribution of Qg is very scattered, especially at the region with large r. For example, for 
a fixed r — 0.22, the values of fig are distributed in a broad range fig E [10^"'^, 10^^"]; 

b. For each fixed r, the values of Qg have an upper limit, and the small region just below this limit tend to 
contain most of the realizations; 

c. For each fixed r, the values of the upper limit Qg from our numerical result are smaller than the analytic 
results of Eq.(25), especially in the region r > 0.01; 

d. At r ~ 0.03, the strength of RGW attains the maximum value fig ~ 2 x 10^^^, which is almost an order 
of magnitude smaller than the analytic result of Eq.(26). And this is beyond the sensitive ranges of LIGO, LISA, 
ASTROD and BBO. 

e. Most of the realizations tend to concentrate in the region with the large values of r, and the larger r is, the 
denser the distribution of the realizations are. More than 90% of realizations are in the region of r > 0.01. This 
phenomenon of distribution may be due to our specific choice of the initial conditions in Eq.(35)-(38). 

Among these 5509 realizations, 50.21% fall into the sensitivity region of the Planck satellite, 97.11% fall into 
that of Clover, 99.29% fafl into that of CMBPol, and 42.91% faU into the sensitivity region of the DECIGO. In 
comparison with the CMB observations, much less realizations are in the sensitivity regions of laser interferometer 
detectors. But DECIGO can detect the RGW with r being much smaller than 10~^, which is beyond the CMB 
experiments. This conclusion is the same as the analytic results in section 3. Therefore, the CMB experiments and 
the laser interferometers are complementary to each other for the RGW detection. 

Now let us look at the 14 realizations that satisfy all the constraints in Eqs.(53)-(56), but end the inflation before 
the e-folds A^ = 70 arrived. We found, for these realizations, the values of e-folds are all in the region of A^ e [40, 70], 
which is consistent with the current observations and the theoretic predictions [39]. In Fig. [3], wc plot them in the 
r — fig plane. This figure shows an interesting feature: larger r corresponds to a smaller fig, which is also consistent 
with the distribution of realizations in Fig. [2]. Among these realizations, 35.71% fall into the sensitive region of 
the Planck satellite, 100% fall into the sensitive regions of Clover and CMBPol, and 64.29% fall into the sensitive 
region of the DECIGO. These results are also consistent with the distribution of the 5509 realizations discussed in 
the above. 

It should be mentioned that, in our numerical calculation, the initial conditions have been chosen randomly for 
the Hubble slow-roll parameters in the regions (35)- (38). It is not clear which one is closer to the actual situation 
of the infiationary process in the early Universe. However, given the very broad range of initial conditions for the 
Hubble slow-roll parameters, these large sample of lO'' realizations may exhaust, to some extent, the reasonable 
reservoir of inflationary models driven by the single scalar field. 

5. Conclusion 

The relic gravitational waves is regarded as a strong evidence for the infiationary models, which is directly related 
to the energy scale of the inflation. Although up to now people have not observed RGW, a lot of constraints 
have been obtained on them. These constraints include two kinds: One is the LSC, which are mainly from the 
CMB observations, especially the recent WMAP results. This can constrain the RGW at very low frequency 
ly G [10~^^, 10~^^]Hz. The other is the SSC, mainly from LIGO, BBN, and pulsar timing, and it is sensitive to the 
waves with high frequencies. A number of experiments arc under development for the RGW detection, which can be 
classified as the following types: the CMB experiments (sensitive to the waves with ly E [10^^^, 10^^''']Hz), including 



Planck, Clover, CMBPol, etc, the laser interferometers (sensitive to the waves with v e [10 ^, 10'']Hz), including 
advanced LIGO, LISA, ASTROD, BBO, DECIGO, etc, the electromagnetic detectors (sensitive to the waves with 
V € [10^ 10i"]Hz). 

We have calculated the strength of RGW, studied how tight of the current constraints on the RGW, and 
examined the detective abilities of the future experiments. When calculating the values of rig(fc), we have worked 
in two methods: the analytic and the numerical. The former method simply shows the dependent relation of ^g{k) 
on the inflationary parameters Ug and r. After considering the current constraints on these parameters, we have 
given an upper limit Q.g < 3.89 x 10~^^, where the three damping effects have been included, such as the reashift- 
suppression, the accelerating expansion, and the free-streaming neutrinos. This limit is in the sensitivity ranges 
of BBO and DECIGO, but beyond which of advanced LIGO, LISA and ASTROD. In the numerical method, the 
energy spectrum ^g{k) has been calculated by solving the inflationary flow equations. The resulting spectrum is 
more precise for the RGW at high frequency range. The corresponding upper limit fig < 8.62 x 10"^'* has been 
given, which is independent of the inflationary parameters and applies to any single- field slow-roll inflationary model. 
After considering the constraints on n^, a and r, this bound becomes fig < 2x 10^^^, which is beyond the sensitivity 
limit of BBO. 

The results from these two methods suggest the consistent conclusions: The current constraints on the RGW 
from LIGO, BBN and pulsar timing are too loose to give any constraint for the single- field inflationary models, 
and the constraints from CMB and LSS are relatively tighter. The future laser interferometer, DECIGO, is more 
effective for detecting the RGW with smaller r, but the CMB experiments, such as Planck, Clover and CMBPol, 
are more effective for detecting waves with larger r. They are complementary to each other for the RGW detection. 
The laser interferometers, as the advanced LIGO, LISA and ASTROD have less chances to find the signal of RGW, 
if the single-field inflationary model is held. 

A final remark should be made, that is, all conclusion on RGW and their detection constraints arrived in this 
paper are pertinent only for single scalar field models for inflation. The RGW generated from other models of 
inflation need to be analyzed separately. 
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Figure 1: The strength of RGW at v = O.lHz depends on the slow-roU parameters Ug and r. This fig- 
ure shows the results of analytic approximation in (25). The solid lines from top down are the curves with 
Us = 1.00,0.98,0.96,0.94,0.90,0.86, respectively. The vertical (dot) lines from right to left are the sensitive limit 
curves of current observations, Planck, Clover, CMBPol, and the sensitive limit of CMB observations, respectively. 
The horizontal (dash) lines from up to down are the sensitive limit curves of BBO and DECIGO, respectively. 
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Figure 2: The distribution of the 5509 reahzations in the r — Vtg plane. The solid lines from top down are the 
analytic curves with n^ = 1.00 and ns = 0.86, respectively. The vertical (dot) lines and the horizontal (dash) lines 
have the same meanings with which in Fig.[l]. 
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Figure 3: The distribution of the 14 realizations in the r — fig plane. The solid lines from top down are the analytic 
curves with n^ — 1.00 and n^ — 0.86, respectively. The vertical (dot) lines from right to left are the sensitive limit 
curves of current observations, Planck, Clover and CMBPol, respectively. The horizontal (dash) lines from up to 
down are the sensitive limit curves of BBO and DECIGO, respectively. 



